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PROCEDE DE SIMULATION BIOMECANI9UE D^TJN ENSEMBLE 
ARTICULATIONS OSSEUSBS 

La present e invention se rapporte au domaine des 
logiciels de simulation biomedicale . 

La pr^sente invention se rapporte plus 
particulierement a un procedg de simulation biom^canique 
d'un ensemble d' articulations osseuses d'un patient, et 
notamment du rachis. Ce proc^dS permet de calculer, estimer 
et visualiser les consequences d'lme operation chirurgicale 
sur une articulation. Le syst^me a ete mis au point pour 
les operations de la colonne vertSbrale et en particulier 
les techniques de stabilisation de la colonne vert^brale. 
D'une manidre g^nSrale, le systdme permet d' informer le 
chirurgien sur I'^tat de I'gquilibre et la repartition des 
efforts dans la colonne du patient avant et apr^s la 
simulation. 



L'art anterieur connait dgjt, par la demande de 
brevet am^ricain US 5 995 738, un dispositif et line m^thode 
pour faciliter 1 ' implantation de composants artificiels 
dans les articulations. L' invention dgcrit des dispositif s 
et des methodes destines H determiner une position 
d' implant pour au moins un composant artificiel dans vine 
articulation et k faciliter son implantation. L' invention 
comprend la creation d'un module de 1 ' articulation du 
patient et la creation d'un moddle du composant a 
implanter. Les modules crees sont utilises pour simuler le 
mouvement de 1 ' articulation du patient selon la position du 
composant. Ce document de l'art anterieur propose done un 
moyen de simulation physique et non virtuel pour analyser 
les mouvements d'une articulation et d'un implant. 

II est propose, dans la demande de brevet americain 
US 6 205 411, un planif icateur de chirurgie assistee par 
ordinateur et un systeme de guidage intra-operatif . 
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L' invention concerne un appareil destine H faciliter 
1' implantation d'un composant artificiel dans une 
articulation. L' appareil comprend un predicteur g^ometrique 
et un simulateur biom^canique du mouvement prS-operatif , 
c'est-a-dire qu'une s^rie de simulations est ef feature sur 
1' implant et 1 ' art iculation avant de proc6der a 
1' operation. 



II est propose, dans la demande de brevet amSricain 
US 5 625 577, une m^thode informatique d' analyse de 
mouvement utilisant la dynamique. L' invention d#crit une 
m^thode permettant d' analyser et d'afficher les mouveraents 
d'un Stre humain. Le corps du sujet est divisS en une 
plurality de segments relics entre eux par des 
articulations. Une fois le corps modSlisg ainsi, il est 
possible de simuler et d' analyser le mouvement du sujet. 

II est ggalement propose, dans la demande de brevet 
PCT WO 99/06960, un systdme et procSdS permettant de 
d^finir et d'utiliser des comportements destines k des 
chaines articul#es dans des animations par ordinateur. En 
1 'occurrence, on difinit pour la chaine articul^e au moins 
une commande telle qu'une commande de forme ou une commande 
de plan de resolution, et on utilise cette commande comme 
une contrainte reprise par le moteur d' animation pour 
1' animation de la chaine avec une cin^matique inverse. 
Chaque commande comporte au moins deux cl^s, chaque cl6 
comprenant un couple constitug d'un vecteur de sens 
effecteur et d'une contrainte associge. Dans le cas des 
clgs de commande de forme, les contraintes associees 
comprennent une liste d' orientations pr6f grentielles des 
membres du corps. Dans le cas des cl6s de commande de plan 
de resolution, les contraintes associees comportent une 
orientation pr§f ^rentielle du plan de resolution. Quels que 
soient les buts assignSs, les cl^s de commande 



wo 2004/061721 PCT/FR2003/003943 

3 

s^lectionn^es siabissent iine interpolation mettant en oeuvre 
des ponderations appropriges visant a I'obtention d'une 
contrainte r^sultante k faire utiliser par le moteur 
d' animation. 

L'art anterieur connait ^galement, par la demands 
de brevet PCT WO 99/60939, un systdme chirurgical 
interactif assiste par ordinateur pour aider le chirurgien 
a positionner des implants dans le femur ou des vis dans 
les pedicules des vertdbres. Ces systSmes apportent une 
assistance ^ la navigation. Cependant, c'est une assistance 
ax#e principalement sur le positionnement , le rep^rage et 
le guidage des ancillaires. Cette technique de navigation 
chirurgicale est assez rgpandue, mais ne permet pas de 
rgpondre k toutes les interrogations du chirurgien. En 
effet, 1' utilisation de ces syst^mes permet de dSfinir la 
trajectoire optimale de la vis p^diculaire, mais ne permet 
pas de savoir par exemple, si elle a 6tS plac6e sur la 
bonne vertdbre. 

Dans la publication scientifique « Morphometric 
evaluations of personalised 3D reconstructions and 
geometric models of the human spine » Aubin C-E et al - 
Medical and biological engineering and computing, Peter 
Perenigrus Ltd. Stevenage, GB - 01/11/1997, plusieurs 
mSthodes de reconstructions ggomgtriques sont gvalu^es. 
Dans la conclusion de cette publication scientifique, les 
auteurs indiquent qu'une des methodes evaluees est 
optimale. L' etude du Dr Aubin confirme qu'une 
reconstruction geom^trique d'une colonne vert^brale a 
partir de plusieurs radiographies constitue une approche 
pertinente et reproductible pour I'^tude et 1' Evaluation 
des dysfonctionnements du rachis ainsi que pour sa 
modelisation par la m^thode des Elements finis. Cette 
publication ne porte aucunement sur les aspects de 
personnalisation mEcanique in vivo (mobility et 
anthropomEtrie) , ni sur la simulation de strategies 
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opSratoires, ni sur leurs methodes d' optimisation par 
approche biomecanique . Enfin I'approche mise en ceuvre dans 
la pr^sente invention associant plusieurs techniques 
(personnalisation geometrique, mobility et anthropometrie 
in vivo, analyse numSrique, critdres d'equilibre, mSthodes 
de recalage entre donnges internes et extemes, calculs des 
efforts) n'y est ni abordee ni sugggr^e. 

Enfin, I'art ant^rieur connalt, par la piiblication 
scientifique « A MRI based semi-automatic modeling system 
for computational biomechanics simulation » Hayasaka T et 
al - Medical imaging and augmented reality, 2001, 
International workshop on 10-12 june 2001 - 10/06/2001, un 
proc^d^ de ggn^ration semi -automati que d'un maillage. Le 
procgdS dScrit dans ce document s' applique au systSme 
cardio-vasculaire, c'est-t-dire aux tissus mous et non aux 
articulations osseuses comme la presente invention. A 
partir d'un module issu d'une base de donn^es et d'une 
image IRM, un algorithme cr6e un module spgcifique en vue 
d'une gtude (nous pr^sumons de la tenue m^canique du tissu) 
bas# sur la mSthode des Pigments finis. L'auteur cite par 
ail leurs des difficultSs restant k r^soudre en vue de crSer 
un maillage utilisable. La simulation m^canique des tissus 
mous qui reste k ce jour du domaine de la recherche 
f ondamentale , est confront^e a la difficult^ de les 
caractgriser m^caniquement in vitro et encore plus in vivo. 
Cette gtude ne porte pas sur la simulation d'un ensemble 
d' articulations osseuses. Elle porte sur un procgd6 de 
generation semi -automati que d'un maillage. 

Ces deux publications scientif iques traitent de 
problgmatiques diffgrentes de celle de la presente 
invention. En particulier, elles ne mentionnent nullement 
d'6tape de personnalisation d'un module numSrique par 
particularisation des paramdtres d' interaction [de 
mobilit^s ou caract^ristiques de rigidit6s] de chacune des 
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articulations reliant lesdits corps rigides en fonction des 
caractSristiques constatSes sur le patient. Cette gtape est 
fondamentale dans la pr^sente invention. 

La prSsente invention entend remSdier aux 
inconv^nients de I'art anterieur en permettant de simuler 
une operation de correction locale ou globale de courbures 
de la colonne vertebrale ou de pose d' implant vertebral a 
partir d' images radiographiques , de series d' acquisitions 
de caractgristiques mesur^es in vivo sur le patient et 
d'une base de donn^es d' implants . La presente invention 
permet de plus de simuler I'etat de la colonne juste aprds 
1' operation chirurgicale . 

Bien entendu, la presente invention ne s' applique 
pas uniquement au rachis. Elle s' applique Sgalement k 
d'autres articulations osseuses comme le genou. 

L' invention consiste en un systSme chirurgical 
assists par ordinateur, permettant au chirurgien de 
simuler, en pr^op^ratoire, les effets sur le patient de la 
chirurgie correctrice qu'il envisage de mettre en 
application. Ce systSme lui permettant de simuler ainsi 
plusieurs strategies opgratoires, il off re au chirurgien un 
outil d'aide au choix de la stratggie opgratoire off rant le 
meilleur compromis entre la stabilisation et la mobility. 

A cet effet, la presente invention concerne dans 
son acception la plus g^nerale un precede de simulation 
biom^canique d'un ensemble d' articulations osseuses d'un 
patient, notamment du rachis, comportant : 

• une gtape d' enregistrement ' d'un module 
num6rique tridimensionnel representee au moins en 
partie par des corps rigides relics par des 
articulations, dans une position de reference ; 
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• une Stape de personnalisat ion de la 
g^ometrie dudit modele [position relative dans 
I'espace de chacvm desdits corps rigides] par des 
donn6es specif iques Sl un patient dans ladite 
position de rgfgrence [par exemple des 
radiographies du patient] ; 

• \ine etape de personnalisation dudit module 
num^rique par part icularisation des paramdtres 
d' interaction [de mobilites ou caract6ristiques de 
rigidity] de chacune des articulations reliant 
lesdits corps rigides en fonction des 
caract^ristiques constat^es sur le patient ; 

caractSris^ en ce que 

I'^tape de particularisation des param^tres 
d' interaction [de mobilites ou caractSristiques de 
rigidity] consiste k : 



• acquSrir les positions dans I'espace d'une partie 
au moins des corps rigides, et a procSder t une 
interpolation pour determiner la position calculSe 
des autres corps rigides pour construire une table 
numgrique comportant les positions relatives de 
chaciin des corps rigides ; 

• exercer au moins une contrainte d^terminee 
sur le patient et a acquerir une information sur la 
position generale d'equilibre r^sultante du patient 
[autre que la position de rgfgrence] ; 

• determiner des fonctions analytiques 
permettant d'approximer les paramdtres 
d' interaction [de mobilit6s ou caractSristiques de 
rigidite] afin de reproduire les positions 
relatives mesurges [fonctions desdites contraintes 
appliquges sur le patient et des modifications 
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geom#triques observ6es sur le patient, et resultant 
de ces contraintes] et ce pour chaque couple de 
corps rigide. 

De preference, le modele num6rique est dgfini par 
des param^tres de positions geom^triques des corps rigides 
et par des paramdtres de rigidite des articulations reliant 
les corps rigides. 

Avantageusement, I'gtape de representation du 
r^sultat d'une contrainte consiste a recalculer le module 
personnalis6 [en position d'gquilibre] resultant d'un 
ensemble de contraintes [par exemple implantation d'une 
prothdse ou d'un implant] comprenant au moins une 
contrainte statique exercSe sur au moins deux corps 
rigides, et imposant un positionnement relatif avec ime 
mobilite ou une raideur diffgrente de celle correspondant k 
la loi de comport ement . 

Selon une variante, I'gtape d' enregistrement du 
module num^rique de 1' ensemble d' articulations standard 
consiste a dgfinir une alternance de corps rigides et 
d' articulations, et k dgfinir pour chacun des couples de 
corps un ensemble de param^tres num^riques caract^risant la 
mobility ou raideur globale resultant de 1' action de 
1' ensemble des glSments intercalaires [par exemple disques 
intervert^braux] et des Elements de liaison [par exemple 
ligaments] ayant un effet sur les parametres d' interaction 
[rigidites] entre les deux corps. 

Selon un mode de mise en ceuvre particulier, l'6tape 
de personnalisation consiste & acquerir au moins une image 
de 1' ensemble d' articulations d'un patient donng, k 
extraire de ladite image les informations nScessaires a la 
construction d'un module r#el par reconnaissance de la 
position des articulations visibles dans ladite image, et 
de modifier le module standard en fonction dudit module 
reel . 
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Avantageusement , I'^tape d' enregistrement d'un 
module numerique consiste a def inir un ensemble standard de 
donnees numeriques comprenant pour chacune des 
articulations representees sous la forme d'un corps 
rigide : 

• un premier descripteur geometrique de 
position de reference, correspondant a la geom^trie 
de 1' ensemble d' articulation pour un patient 
« standard » dans une position <cde rgf6rence», 
ledit descripteur 6tant determine pour chaque corps 
rigide de fagon relative par rapport k un corps 
ad j acent ; 

• un deuxieme descripteur m^canique 
d' interaction avec chacun des corps adjacents, 
ledit descripteur m^canique 6tant representatif de 
la loi de comportement lorsque 1' ensemble 
d' articulation est soumis a 1' action d'au moins une 
contrainte exterieure ; 

I'^tape de personnalisation consistant a modifier 
ledit ensemble standard de donnees par des donnees 
personnalisSes . 

Selon un mode de mise en oeuvre particulier, le 
procgdg comporte en outre une etape de recalage consistant 
a mettre en correspondance des donnees d' images radios et 
des donnees d ' acqui sit ion externe, cette 6tape se 
decomposant en deux sous -Stapes : 

• recalage de la reconstruction radio par 
rapport a la courbe 3D issue de donnees 
d' acquisition externe dans la m§me 
position ; 

• determination de la repartition des points 
de la courbe 3D associ^s aux vertdbres. 
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positionnSs dans le repdre de Stokes et leur 
tangente associ6e. 

On comprendra mieux 1' invention a I'aide de la 
description, faite ci-apres a titre purement explicatif, 
d'un mode de realisation de 1' invention, en r^fSrence a\ix 
figures annexees : 

• la figure 1 illustre 1 ' architecture de la 
personnalisation mecanique du module ; 

• la figure 2 repr^sente 1' architecture g^nSrale du 
simulateur selon un mode de realisation de 
1 ' invention ; 

• la figure 3 reprgsente 1 ' architecture du module 
selon un mode de mise en oeuvre de 1' invention ; 

• la figure 4 illustre la mise en oeuvre du proc^d^ 
de recalage ; 

• la figure 5 illustre la definition des angles 
intervertebraux ; 

• la figure 6 repr^sente un exemple des relations 
entre les angles intervert^braux dans le plan 
frontal (abscisse en •> et ordonnSe en *>) ; 

• la figure 7 illustre le calcul des centres de 
rotation ; 

• la figure 8 represent e un exemple de calcul des 
centres de rotation sur un sujet scoliotique ; 

• la figure 9 illustre 1 ' interpolation des 
couronnes ; et 

• la figure 10 repr^sente \m modele de repartition 
des masses dans le tronc. 

L' invention se distingue des produits connus de 
guidage chirurgicaux, utilises pour assister le chirurgien 
pendant son operation. En effet, 1' invention consiste 
notamment a modeiiser la colonne vertebrale de I'individu 
qui va gtre opere, de simuler la pose de 1' implant ou de 
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prothdse et de calculer la position d'6quilibre de 
I'individu une fois 1' implant ou la prothdse pose. 

Le traitement des donnees etant logiciel, 
1' invention concerne aussi 1 ' architecture logicielle mise 
5 en place pour la realisation de 1 ' architecture 
f onctionnelle. La structure logicielle comprend plusieurs 
serveurs de bases de donnees : une base de donnees 
utilisateurs, une base de donnees vertdbres, une base de 
donnees patients et une base de donnees implants. Les 

10 donnees contenues dans ces bases sont interrog^es et mises 
a jour par differents utilisateurs afin de construire le 
modele 3D de la colonne du patient puis de simuler les 
consequences de la mise en place d'un implant. 

Le systdme selon 1' invention est un simulateur des 

15 consequences biom^caniques et cinSmatiques de 1 ' acte 
chirurgical des traitements des pathologies rachidiennes 
sur la morphologie du patient. La Figure 2 represente 
1 ' architecture ggngrale du simulateur. 

II doit permettre k un chirurgien d'optimiser et 

20 ameliorer la planif ication de son traitement. Ce simulateur 
pourra lui apporter une meilleure connaissance des 
propriStes geomStriques et m^caniques des divers tissus de 
la colonne vert^brale. II pourra aussi tester diff^rentes 
approches de son geste pour permettre une correction 

25 optimale. 

II doit rSpondre k un besoin de sante, puisque la 
tendance est a la recherche du confort, de la securite, de 
la quality et de la fiabilite. L'alldgement des traitements 
mgdicaux se traduit par une meilleure quality de vie post- 
30 opgratoire et en particulier par des hospitalisations de 
frequence moindre et de duree plus courte. 

Enfin, il peut egalement avoir des retombees dans 
le domaine de 1 ' enseignement pour 1 ' apprentissage des 
futurs chirurgiens. 

35 
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Le simulateur est un outil d'aide a 1' analyse de la 
faisabilitS du geste. II permet de simuler I'Squilibre pre-- 
operatoire, les efforts inter-segmentaires correspondants 
tenant compte de I'effet des rigidites musculaires et 
ligamentaires, apprehender 1' evolution de cet equilibre et 
de ces efforts en post-operatoire en fonction des courbures 
introduites par le chirurgien lors de 1' operation. Le 
chirurgien devra indiquer en entree la deform^e voulue de 
la zone instrunientee . Pour cela une approche dite 
« globale » sera utilisSe. 

Deux radiographies sont n^cessaires pour la 
realisation d'une simulation. Les images a a j outer au 
systdme peuvent Stre : 

- le fichier rSsultat de la num6risation (scan) 
d'une radiographic traditionnelle ; 

- un fichier fourni par un autre appareil de 
radiologie ( radiographic numerique) . 

L' alimentation du systdme peut §tre manuelle (un 
utilisateur va fournir les fichiers d' images au systeme) ou 
automatique, les images sont stock^es directement par les 
appareils de radiologie et recup^rges via un reseau 
intranet/ internet . 

Deux sortes de numSrisations seront possibles : 

- numgrisation manuelle : des points specifiques 
doivent Stre reper^s manuellement par un acteur sur ces 
images ; 

- numSrisation semi-automatique . 

Les points specifiques detectes par la numerisation 
vont permettre de calculer les coordonnees 3D des 
vertebres (donnees geomStriques) . 

Les donnees geometriques issues des images 
numeriques vont permettre de construire un modele de la 
colonne du patient en trois dimensions, Ce modSle rSsulte 
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de 1' adaptation d'un modele 3D standard, predefini dans le 
systdme, aux caracteristiques ggometriques du patient. 

L'utilisateur doit pouvoir visualiser le modele 3D 
du rachis selon les plans frontal, sagittal et apical ; de 
plus, il doit pouvoir le comparer aux radiographies 
utilis^es pour sa construction, 

Les caracteristiques m^caniques du patient 
(resultats issus d' acquisitions (tests cliniques) ) 
permettent de personnaliser le moddle. Pour cette 
operation, un modele g^ometrique standard de la colonne 
vertebrale est enrichi par les donn^es m^caniques 
personnelles du patient conceme. 

Les caracteristiques biomecaniques de la colonne du 
patient (angle de scoliose, rotation axiale, etc) sont 
calculSes k partir des radiographies du patient (avant la 
premiere simulation), ou k partir du modSle rSsultant d'\me 
simulation. Certains paramStres pourront Stre pr^-calcul^s 
(pente sacrSe, angle d' incidence, gite sagittale, courbures 
rachidiennes) sur les radiographies. 

Lors de la simulation, l'utilisateur doit avoir la 
possibility de choisir des segments de colonne sur lesquels 
il impose des d^placements . Lors de la simulation de 
I'gquilibre, l'utilisateur va pouvoir entrer les valeurs de 
certains paramdtres cliniques qu'il souhaite simuler pour 
un segment donnS. L'utilisateur doit pouvoir consulter les 
mobilitSs ou rigidites du rachis, la graduation sera 
normee, Le type de la tige utilisee pour deformer un 
segment peut etre choisi parmi plusieurs propositions 
(rigidity, diamdtre, etc.)* 

A la suite d'une manipulation du moddle par 
l'utilisateur, le logiciel doit verifier la validite des 
actions effectuees, et prevenir l'utilisateur de toute 
incoherence (valeur impossible a obtenir) . 

Le simulateur doit afficher les nouvelles courbures 
du modele et la nouvelle position (stature) de l'6quilibre 
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du patient. II doit Stre possible de comparer la (les) 
courbure(s) avec la (les) courbure(s) initiale(s) (avec le 
module 3D ou la radiographie) . 

L' Evolution relative des efforts intervert^braux 
5 est quantifige et representee graph i quement . 

L'utilisateur doit pouvoir comparer les rigidit^s 
du rachis avant et aprds la simulation. 

Les efforts resultants dans la tige suite aux 
manipulations effectuees doivent indiquer si la tige va se 
10 d^former. 

Des dossiers patients, contenant les informations 
concemant la simulation, doivent pouvoir atre cr6Ss par un 
utilisateur et maintenus §l jour par le systdme, (sauvegarde 
de plusieurs simulations dans xin mSme dossier) . De plus, le 
15 systdme doit pouvoir retirer ou saisir les informations 
concemant un patient (Sge, poids, taille, etc.) dans des 
systdmes de donriSes externes, sous reserve de leur 
existence. 

Les droits et profils utilisateurs vont permettre 
de diff^rencier les fonctionnalitgs disponibles pour chaque 
utilisateur. En fonction de leurs droits et profils, les 
utilisateurs auront accds $l des interfaces graphiques 
diff#rentes. De plus, un systdme de prgfgrence pourra §tre 
mis en place pour chaque utilisateur (position des menus, 
25 page d'ouverture etc.), ainsi qu'un mecanisme de 
traqabilite, qui permettra de suivre les evolutions d'un 
dossier patient. 

A tout moment, et en fonction de ses droits, un 
utilisateur peut avoir accSs aux donnSes d' entries, c'est- 
a-dire visualiser, remplacer ou modifier les images 
utilisees pour la construction du module, les points 
numeris^s ou les operations effectuees pour le traitement 
semi- automat ique de ces images, les implants places sur le 
moddle, les caract6ristiques m^caniques du patient... 



20 



30 
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L'historique d'une simulation represente la 
totalite des actions effectuees sur un patient (paramdtres 
de la simulation) . L' utilisateur doit pouvoir revenir sur 
une action passee (changer une radiographie, ressaisir les 
points pour la numerisation, calculer de nouveaux 
parametres cliniques, ou modifier les actions faites au 
cours de la simulation) 

La Figure 1 represente 1 ' architecture de la 
personnalisation mecanique du moddle. 

La personnalisation mecanique du modSle est bas6e 
sur trois types de donnees : 

- les radiographies du patient avec des repdres 
cutanes ; 

- 1' acquisition de la courbure generale du rachis 
sous diffgrentes postures caracteristiques ; 

- des donnees anthropom^triques . 

A partir de ces donnges, des traitements utilisant 
les lois de la mecanique nous permettent d'obtenir : 

- la geometrie du patient ; 

- les paramdtres cliniques ; 

- le moddle mecanique personnalisS . 

La numerisation des radiographies du patient nous 
permet d'obtenir une precision de six points par vertdbre. 
En effectuant une extrapolation, on arrive k une precision 
de douze points par vertebre , 

L'etape suivante est 1 ' acquisition de la courbure 
generale du rachis sous diffgrentes postures 
caracteristiques • 

Pour ce faire, lors d'un examen clinique, le 
patient etudig est soumis a une serie de tests au cours 
desquels sera evaluSe la ligne genSrale de la colonne au 
niveau dorsal. Le positionnement du bassin et des 6paules 
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sera n€cessaire pour dgfinir les orientations des 
extr^mit^s de la colonne. 

Le patient devra Stre, eventuellement , maintenu au 
niveau du bassin pour limiter 1 ' intervention des membres 
exterieurs dans 1' gtablissement de son equilibre general, 
il sera effectu^ une acquisition k I'aide d'un appareil, 
des positions dans I'espace de repSres cutan€s 
identif iables lies aiix vertdbres scoliotiques par rapport a 
^xn. r^f^rentiel connu. 

Une acquisition sera effectuge patient au repos 
dans des conditions aussi proches que possible des 
conditions de prise de la radiographie calibr^e patient 
debout au repos. Cette acquisition sera utilis^e pour 
dSterminer la ligne des centres des corps vert^braiix a 
partir de la ligne ggnSrale de la colonne. Ceci est rendu 
possible par le fait que la radio calibrge peut donner les 
coordonn^es de tous ces points gr^ce aux repdres plombSs . 
Le positionnement des vertdbres peut done Stre d^terming en 
fonction des repSres cutan^s pour une position donnge, 
debout, au repos. II sera ensuite gtabli \me transformation 
corrigge entre les repdres cutan^s et le positionnement des 
vert^bres prenant en compte 1' influence de la cin#matique 
des vertabres entre elles lors des mouvemehts du patient. 

A partir des donn^es anthropom^triques , on 
determine les masses segmentaires et les centres de masses, 
permettant ensuite de determiner les moments et les centres 
de masse par niveau vertebral . 

L' action musculaire sera dissoci^e de 1' action 
inter segmentaire. 



Les donn^es disponibles dans la litt^rature 
permettent d'4valuer qualitativement la forme des lois de 
comportement admissibles pour le module. 
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Les lois de comportement doivent repondre aiix 
exigences suivantes : 

la plupart des lois doivent suivre des 
comport ements impairs ; 

- le comportement asymptotique doit St re assure ; 

- les ph^nomenes de couplage sont pris en comptes ; 

- la pertinence et la simplicity de calcul. 

Les lois de comportement devront Stre recalculges 
pour chaque niveau vertebral et pour chaque patient afin de 
prendre en compte les singularit^s dues aux pathologies 
6tudi4es . 



Selon une variante de 1' invention, on met en 
correspondance les radios avec les acquisitions extemes. 
La mise en correspondance se passe en deiix Stapes : 

• Le recalage de la reconstruction radio par rapport a 
la courbe 3D issue de donnes d' acquisition exteme 
dans la m§me position. 

• La determination de la repartition des points de la 
courbe 3D associSs aux vertdbres, positionnes dans le 
repSre de Stockes et leur tangente associees. 

Lors des acquisitions radio et de donnees externes, trois 
repSres spatiaux sont Sgalement enregistres. Les deux 
extrSmitSs de la rSglette, et une bille en plomb (point 
bas) , plac6e en bas de la colonne. 



Ces trois points dont on connait la position 
relative, k la fois dans le repere de reconstruction radio 
et dans le repere des donnSes d' acquisition extemes, 
permettent de recaler les deiix enregistrements dans le m§me 
repere, celui de la radio. 
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Les trois points sont suffisants a determiner 
un repere commun dans lequel 1' ensemble des donnees est 
connu, ce qui permet done de definir la matrice de passage 
entre les deux enregistrements . 

Compte tenu des erreurs d' acquisition possible, 
rotation du poignet involontaire, rugosite de la peau 
variable le long de 1 ' enregistrement , trois autres 
recalages sont egalement effectues : 

1) on compare la distance entre le point 
bas de la r^glette et le point bas plac6 sur la peau, 
pour la radio et 1 ' acquisition externe. La difference 
entre les deux enregistrements quant if ie le trajet 
effectug par la spline sur la peau et done rea juste 
1' acquisition vis-S-vis de la pression exercee par 

15 I'operateur sur la peau du sujet. 

2) Cet ajustement est complStS par un 
ajustement des longueurs des acquisitions externes de 
telle maniSre k garantir le fait que les quatre 
acquisitions partent du mSme point. 

3 ) Lors de 1 ' acquisition, suivant le 
positionnement de I'op^rateur, il est frequent de 
constater un biaisement de 1 ' acquisition vers la 
gauche ou la droite. Ce decalage est corrige en 
redressant les acquisitions en contraignant 
1' acquisition allant jusqu'au point bas de respecter 
la jonction entre le point concourant defini au 2 et 
le point bas defini sur la radio. 

Suite a ces trois recalages, il est possible 
d' avoir une confiance certaine dans le recalage effectue, 
30 corrige des erreurs. 

Repartition des vertdbres sur la courbe 3D : 
A partir du recalage entre les radios, la 
deuxidme 6tape consiste S placer les vertdbres sur la 
35 spline. On considdre que le point de correspondance entre 
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la vertSbre et la spline est le point d' intersection du 
repere de Stockes (plan XY) avec celle-ci, 

I>our les vertebras du bas de la colonne, la 
distance importante entre la colonne et la surface du dos, 
ainsi que la forte inclinaison des vertdbres provoquent des 
posit ionnements non successifs des vertdbres sur la Courbe 
3D. Sur cette portion d' enregistrement , on applique la 
demarche inverse, on place des points gquidistants puis on 
definit leurs positions dans le repdre de Stokes de la 
vertSbre correspondante . 

On enregistre ^galement les tangentes de la 
Courbe 3D auix points consideres de maniSre t pouvoir par ia 
suite repositionner le modele. 
15 Cette mgthode permet de dSfinir egalement sur 

les enregistrements d' acquisitions externes la position du 
sacrum et 1' angle du bassin avec la verticale, appele angle 
sacr^ , 

Meme si ces parties anatomiques ne font pas 
parties de la reconstruction k partir de la radio, ils sont. 
utiles pour dSfinir I'Squilibre du rachis. 

On suppose ensuite que la repartition des 
vertdbres sur toutes les acquisitions externes reste la 
m§me sur toutes les acquisitions, ce qui permet de 
positionner les vertdbres dans toutes les positions. 



20 



25 



Notons quelques remarques sur 1' analyse des donnSes : 

o L' etude se fait sur les deux plans 
frontal et sagittal separSment . En effet, le choix a 
ete pris d'Studier dans un premier temps la colonne 
vertebrale de manidre dScouplee. C'est pourquoi les 
mouvements demand^s aux patients sont des inflexions 
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laterales pour 1' analyse du plan frontal et des 
flexions/extensions pour 1' analyse du plan sagittal. 

Le calcul des angles intervert6braux est 
5 illustre sur la figure 5. 



Le but est de trouver une relation liant les 
angles intervertebraux. Pour cela, un module d' analyse a 
6te ajoutg, tragant un angle intervertebral en fonction de 

10 celui qui se trouve « en dessous » lorsque I'on parcourt la 
colonne de bas en haut. 

Exemple de trac6 des relations entre les angles 
intervertebraux (le premier graphique reprSsente 0(C7.n)/(nr2) 
haut k gauche, puis on parcourt la colonne en descendant 

15 idans le sens de lecture classique) :La figure 6 illustre un 
exemple des relations entre les angles intervertebraux dans 
le plan frontal (abscisse en ^ ; ordonnee en ^) 

Les conclusions de ces nuages de points sont 
20 les memes pour tous les sujets. La distribution est 
consideree lineaire. Le modele est done batit sur ces 
constatations : 

On definit 1' Equation de regression sur les 
points a chaque niveau intervertebral. A chaque equation on 
adjoint des bornes caract^risant la nature bornee de la 
mobilite d'une vertebre par rapport a celle qui se trouve 
en dessous. On fait alors I'hypothdse que les mouvements 
d' inflexion laterale et de flexion/extensions maximales 
definissent ces butSes. 



25 



30 



35 



Cette approche a 6t6 validee en analysant un 
certain nombre de sujets sains. Pour la plupart des sujets, 
et des niveaux intervertebraux, la modeiisation adoptee 
semble tout a fait appropriee : les coefficients de 
correlation sont le plus souvent superieurs a 0,8. 
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Cependant, il existe des niveaux pour lesquels la 
modSlisation est incorrecte. Ce phenomene correspond k des 
niveaux bloques ou presque bloqu^s chez le patient. On 
considSre done qu'il reste correct de construire la 
regression lineaire, et les but^es contraindront de ce fait 
suf f isamment le moddle. 

Cette methode permet de definir des lois de 
comportement du disque T1/T2 a L5/S1. 

Les centres de rotations d'une vert^bre par 
rapport ^ une autre se calculent k partir des diffgrentes 
positions enregistrges dans le plan frontal (bending) et 
sagittal (Flexion/extension) . Les points de la vertdbre 
supgrieure correspondant ^ chaque position sont places dans 
le repSre de Stockes de la vertdbre inf^rieure. Ces points 
forment alors une trajectoire consid^r^e comme circulaire 
dont on calcule le centre de rotation par la mSthode des 
moindres carr^s. 

Ceci est illustrS sur la figure 7. 

Cette m6thode permet d' identifier des centres 
de rotation apparents et personnalis^s verifiant les 
donn€es exp^riraentales obtenues sur un patient. 

Un exemple de calcul des centres de rotation 
sur un sujet scoliotique est illustre sur la figure 8. 

On constate que cette mgthode permet de 
retrouver la lordose lombaire Sl partir des acquisitions 
externes. Cette methode n^cessite un minimum de trois 
acquisitions distinctes . 

Ces centres de rotations sont calcules dans le 
rep^re global de la radio. Pour chaque vertebre, deux 
centres de rotation distincts sont calculus dans I'espace, 
un pour les mouvements sagittaux et un pour les mouvements 
frontaux. Ces centres sont ensuite exprim^s dans le repdre 
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de Stokes de la vertebre infer ieure de manidre a assurer 1 
repositionnement dans le solver. 



Un exemple d'algorithme est represente sur 1 
page suivante . 



wo 2004/061721 



22 



PCT/FR2003/003943 



• Acquisition du module cindmatique : 

ModMe gdometrique issue des radios 

Relation lineaire entre les differents angles intervertebraux 

Amplitude maximum des diffdrents angles intervertebraux 

Position dans les rcpferes de Stokes des vertfebres des centres de rotation frontaux et sagitta 



• Initialisation du variant : angle intervertebral L5/S 1 

• Calcul des angles intervertebraux initiaux sur le modele g^om^trique 

Tant que le crit^re de confort n'est pas optimise 
Faire 

L' angle variant est incremente du pas de calcul (par defaut 0,001) 

Si Tangle intervertebral variant depasse ses butees Alors Tangle vaut la valeur de la butee 
et le variant est Tangle suivant 



Pour tons les angles suivants 
Faire 

L' angle intervertebral i est incremente en fonction de Tangle i+1 

Si Tangle intervertebral i depasse ses butees Alors Tangle vaut la valeur de la butee 
Fin Faire 

Repositionnement du module a partir des nouveaux angles intervertebraux. 
Application du critere mecanique pour repositionnement du module autour de Taxe des 
t8tes femorales 

Recalcule des angles intervertebraux sur le modfele geometrique suite i sa deformation 
Calcul du critere de confort 

Si le niveau de confort dinainue, on incremente le modele dans T autre sens. 

Si les deux sens de parcours sont faits, on diminue le niveau d' exigence du critere de 

confort 



Fin Faire 
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La colonne vertSbrale est jusqu'^ maintenant 
simulSe cin#matiquement . On sdDoutit k une g^omStrie du 
rachis repr^sentant I'gquilibre du patient. II est 
5 int^ressant d^sormais de connaitre la repartition des 
efforts qui s'exercent sur la colonne, afin de construire 
un moddle m^canique complet. L'art ant^rieur connalt une 
approche pour le calcul des efforts s'inspirant du modele 
anthropom^trique . Le principe de cette approche est de 

10 d^couper le tronc en quatre tranches rattachees a quatre 
parties bien precises de la colonne vertebrale et de 
calculer leur poids et centre de gravite. Le contour de 
chaque tranche est obtenu par dimensionnement de patron 
g^n^rique de la population moyenne aux dimensions du sujet 

15 mesure . 

La m^thode utilisee dans le moddle selon 
1' invention est sensiblement diffgrente. Cinq couronnes 
sont enregistrges a partir d' enregistrements de contours. 
Ces couronnes sont ensuite discr^tisees en 60 points. Dans 
20 un deuxidme temps, ces couronnes sont positionnges par 
rapport a la colonne k partir du recalage. Ce 
positionnement se fait en confondant le point sur la peau, 
au niveau de la vertebre associ^e k la couronne, enregistrg 
par la courbe de recalage avec le point de la couronne 
correspondant au centre du dos. Les autres couronnes sont 
ensuite inteipolSes pour obtenir 18 couronnes d^finissant 
17 tranches centr^es verticalement sur le centre du corps 
vertebral de la vertdbre correspondante . 

La Figure 9 illustre 1' interpolation des couronnes. 
Le volume et le centre de gravity de chaque tranche ainsi 
d^finie est calculi en les d6coupant en prismes 
Sl€mentaires . 
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A partir des volumes et des centres de gravit^s de chaque 
tranche, on dresse le module suivant : 

5 Le poids de la tSte et des bras est repris uniquement 

par la colonne vertebrale 

Chaque tranche est divis^e en deux : les viscSres, et 
les parties solides comme les cotes, muscles, peau 
etc . 

10 o La fraction solide de la tranche est reprise 

entiSrement par la colonne 
o La fraction visc^re a un comportement 
hydrostatique : 

■ Les viscSres n'engendrent pas de moment 

■ Seule la composante normale a la colonne 
vertebrale voit les efforts transmis par 
les viscdres 

■ Les composantes tangentielles k la colonne 
sont transmises k la tranche d'au-dessous 
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La Figure 10 reprgsente un moddle de 
repartition des masses dans le tronc. 

On obtient alors le torseur des efforts sur chaque vertSbre 
dependant uniquement de la g^om^trie entree. 
La repartition visc^re / corps dur pour chaque tranche est 
pour 1' instant defini a partir d'un nombre restreint de 
donnees. Cette repartition sera par la suite affinee. 

L' architecture du programme se base sur une 
structure modulaire correspondant Sl chaque etape de calcul 
independant, de mani^re H pouvoir faire evoluer le 
programme sans alterer son f onctionnement . Les modules de 
calcul sont prgsentes ci-dessous. lis sont classes en 4 
categories : 

• Les modules d' acquisitions 

• Les modules de definition du module cinematique 
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• Les modules du solveur de I'Squilibre 

• Module de calcul des efforts sur le rachis 
instrximente . 

5 L' invention est decrite dans ce qui pr^cdde ^ titre 

d'exemple. II est entendu que I'homme du metier est a mSme 
de rSaliser diffSrentes variantes de 1' invention sans pour 
autant sortir du cadre du brevet. 
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1 - Proc^dg de simulation biomecanique d'un 
ensemble d' articulations osseuses d'un patient, notamment 
5 du rachis, comportant : 

• une gtape d' enregistrement d'un module 
numerique tridimensionnel representee au moins en 
partie par des corps rigides relics par des 
articulations, dans une position de reference ; 

• etape de personnalisat ion de la 
g^ometrie dudit module [position relative dans 
I'espace de chacun desdits corps rigides] par des 
donnges spScifiques a un patient dans ladite 
position de reference [par exemple des 

15 radiographies du patient] ; 

• une etape de personnalisation dudit moddle 
numerique par particularisation des paramdtres 
d' interaction [de mobilites ou caracteristiques de 
rigidite] de chacune des articulations reliant 

^° lesdits corps rigides en fonction des 

caracteristiques constatees sur le patient ; 

caracterise en ce que 

1' etape de particularisation des parametres 
d' interaction [de mobilites ou caracteristiques de 
rigidite] consiste H : 



- acquerir les positions dans I'espace d'une 
partie au moins des corps rigides, et ^ proceder a 
une interpolation pour determiner la position 
calcuiee des autres corps rigides pour construire 
une table numerique comportant les positions 
relatives de chacvm des corps rigides ; 
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• exercer au moins une contrainte d#terminge 
sur le patient et a acquerir une information sur la 
position generale d'equilibre r^sultante du patient 
[autre que la position de rSf^rence] ; 

• determiner des fonctions analytiques 
permettant d'approximer les paramStres 
d' interaction [de mobilit^s ou caractSristiques de 
rigidite] afin de reproduire les positions 
relatives mesur^es [fonctions desdites contraintes 
appliquees sur le patient et des modifications 
geom^triques observees sur le patient, et resultant 
de ces contraintes] et ce pour chaque couple de 
corps rigide. 



2 - Proc#dS de simulation biom^canique d'un 
ensemble d' articulations osseuses selon la revendication 1, 
caract^risg en ce que le module num^rique est dSfini par 
des paramStres de positions ggomStriques des corps rigides 
et par des paramStres de rigiditg des articulations reliant 
les corps rigides. 

3 - Proc^de de simulation biomecanique d'un 
ensemble d' articulations osseuses selon la revendication 1 
ou 2, caractSrisg en ce que I'^tape de representation du 
rSsultat d'une contrainte consiste a recalculer le module 
personnalisg [en position d'gquilibre] resultant d'un 
ensemble de contraintes [par exemple implantation d'une 
prothese ou d'un implant] comprenant au moins une 
contrainte statique exercge sur au moins deux corps 
rigides, et imposant un positionnement relatif avec une 
mobility ou une raideur diffgrente de celle correspondant a 
la loi de comportement . 
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4 - Proc^dS de simulation biomecanique d'un 
ensemble d' articulations osseuses selon les revendi cat ions 
1 et 2, caract#risS en ce que l'6tape d' enregistrement du 
module numgrique de 1' ensemble d' articulations standard 
consiste k dgfinir une alternance de corps rigides et 
d' articulations, et a d^finir pour chacun des couples de 
corps un ensemble de paramdtres numSriques caracterisant la 
mobility ou raideur globale resultant de 1' action de 
1' ensemble des Pigments intercalaires [par exemple disques 
intervert^braux] et des Pigments de liaison [par exemple 
ligaments] ayant un effet sur les paramdtres d' interaction 
[rigidit^s] entre les dexix corps. 

5 - Proc^dS de simulation biomScanique d'un 
ensemble d' articulations osseuses selon I'une au moins des 
revendications prec^dentes, caract^risg en ce que I'^tape 
de personnalisation consiste k acqu^rir au moins une image 
de 1' ensemble d' articulations d'un patient donn6, ^ 
extraire de ladite image les informations n^cessaires k la 
construction d'un module r€el par reconnaissance de la 
position des articulations visibles dans ladite image, et 
de modifier le moddle standard en fonction dudit module 
r€el . 



6 - Proc^dg de simulation biomecanique d'un 
ensemble d' articulations osseuses selon I'vme au moins des 
revendications prScedentes, caract^rise en ce que I'gtape 
d' enregistrement d'un modele numerique consiste ^ dSfinir 
un ensemble standard de donnees numeriques comprenant pour 
chacune des articulations representees sous la forme d'un 
corps rigide : 

• un premier descripteur g^om^trique de 
position de rgfgrence, correspondant k la ggomgtrie 
de 1' ensemble d' articulation pour un patient 
« standard » dans une position «de reference*. 
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ledit descripteur gtant determine pour chaque corps 
rigide de fagon relative par rapport k un corps 
ad j acent / 

• un deuxidme descripteur m^canique 
d' interaction avec chacun des corps adjacents, 
ledit descripteur mecanique 6tant repr^sentatif de 
la loi de comportement lorsque 1' ensemble 
d' articulation est soumis a 1' action d'au moins una 
contrainte exterieure ; 
I'^tape de personnalisation consistant a modifier ledit 
ensemble standard de donn^es par des donn^es 
personnalis^es . 



7 - Proc^dg de simulation biomScanique d'un 
15 ensemble d' articulations osseuses selon I'une au moins des 
revendications pr6c6dentes, caractSrise en ce qu'il 
comporte en outre une ^tape de recalage consistant k mettre 
en correspondance des donnges d' images radios et des 
donnges d' acquisition exteme, cette 6tape se d^composant 
20 en deux sous -Stapes : 

• recalage de la reconstruction radio par 
rapport k la. courbe 3D issue de donnees 
d' acquisition externe dans la meme 
position ; 

• determination de la repartition des points 
de la courbe 3D associes aux vertdbres, 
positionn^s dans le repere de Stokes et leur 
tangente associSe . 
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Figure 1 

Architecture de la personnalisation mecanique du modele 
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Figure 2 
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Figure 5 
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Figure 6 
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Figure 7 
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Figure 9 
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Figure 10 
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